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C-Gonsrdter  the  Wmimar  tin  gale  regression  model 


X(t)  =  A'(0)+  f  Y{9)a{3,Z{3))ds  +  M{t), 
Jo  . . 


where  Y,  Z  are  observable  covariate  processes,  ^  is  a  (deterministic)  function  of  both,  time  and 
the  covariate  process  Z,  and  Af  is  a  square  integrable  martingale.  Under  the  assumption  that  i.i.d. 
copies  of  X,  Y,  Z  are  observed  continuously  over  a  finite  time  interval,  inference  for  the  function 

A  1*  ^ 

-n{t,z)  is  investigated.  An  estimator  A  for  the  time  integrated  A{t,z)  and  a  kernel  estimator  of 
a(t,  z)  itself  are  introduced.  For  X  a  counting  process,  A  reduces  to  the  Nelson-Aalen  estimator 
when  Z  is  not  present  in  the  model.  Various  forms  of  consistency  are  proved,  rates  of  convergence 
and  asymptotic  distributions  of  the  estimators  are  derived.  Asymptotic  confidence  bands  for  the 
time  integrated  and  a  Kolmogorov-Smirnov-type  test  of  equality  of  a  at  different  levels  of 

the  covariate  are  given. 
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1.  Introduction. 


A  useful  way  of  modelling  the  dependence  of  a  counting  process  iV  on  a  covariate  process  Z 
was  given  by  Aalen  (1975).  In  his  multiplicative  intensity  model  N  is  supposed  to  have  an  intensity 
\  given  by 


A{t)  =  a(t)  Z{t), 

where  a  is  an  unknown,  deterministic  function  of  time  (the  hazard  function).  In  the  present  paper 
we  study  inference  for  counting  processes  with  intensities  having  general  dependence  on  a  covariate 
process  Z,  as  in 


X{t)  =  Y{t)a{t,Z{t)),  (1.1) 

where  a  is  an  unknown,  deterministic  function  of  both  time  and  the  covariate  process  Z.  The 
covariate  process  Y  is  taken  to  be  an  indicator  process,  assuming  the  value  1  when  the  counting 
process  is  under  observation,  zero  otherwise. 

An  important  example  of  our  model  arises  in  survival  analysis.  Suppose  that  the  conditional 
hazard  function  h{t\Z)  for  the  survival  time  T  of  an  individual  given  the  covariate  process  Z 
has  the  form  h(t\Z)  =  a{t,Z{t)).  The  observable  portion  of  the  individual’s  lifetime  is  given  by 
T  =  min(T,  C),  where  C  is  a  censoring  time.  We  observe  T,  6  =  I[T  <  C)  and  Z[t)  for  t  <T.  Let 
N{t)  =  /(T  <t,6  —  1),  the  counting  process  with  a  single  jump  at  an  uncensored  survival  time.  If 
T  and  C  are  conditionally  independent  given  Z  then  N  has  intensity  (1.1),  where  Y[t)  =  /(T  >  t) 
is  the  indicator  that  the  individual  is  “at  risk”  at  time  t.  We  shall  introduce  an  estimator  A  for 
the  time-integrated  conditional  hazard  function  A(*,  z)  =  a(s,  z)  ds  which,  in  the  special  case 
of  a  time-independent  covariate  Z,  coincides  with  an  estimator  proposed  in  unpublished  work  of 
Beran  (1981).  Dabrowska  (1987)  recently  obtained  a  weak  convergence  result  for  Beran’s  estimator 
by  proving  a  “conditional”  ansJogue  of  Theorem  4  of  Breslow  and  Crowley  (1974).  We  obtain 
asymptotic  results  for  our  estimator  by  using  a  martingale  approach  (in  particular,  Rebolledo’s 
martingale  central  limit  theorem)  which  enables  us  to  give  quite  simple  proofs  and  to  avoid  the 
restrictive  assumption  of  time-independent  covariate  Z. 

For  another  example  of  our  model,  consider  a  pure  jump  process  describing  the  motion  of  a 
particle  on  a  finite  state  space  (1,2,  ...,m}.  Let  the  intensity  Qjy(t,s)  of  trauisition  from  state  »  to 
state  j  depend  on  the  (calendar)  time  t  and  on  the  time  s  spent  in  state  t  since  the  leist  jump.  Let 
yi(t)  be  the  indicator  that  the  particle  is  in  state  i  at  time  t.  Then  the  counting  process  iV,y(t) 
which  registers  the  number  of  transitions  from  state  i  to  state  j  up  to  time  t  has  intensity 

X{t)  =  Yi{t-)ai,{t,L{t-)),  (1.2) 

where  L{t)  is  the  length  of  time  which  at  time  t  has  elapsed  since  the  last  jump.  In  the  terminology 
of  Markov  renewal  processes  (Pyke,1961),  L{i)  is  the  backward  recurrence  time.  In  the  ceise  that 


a,-y  only  depends  on  calendar  time  t,  inference  for  onj  has  been  studied  by  Aalen  (1975,  1978).  In 
the  case  that  oiij  only  depends  on  the  backward  recurrence  time  L{t),  Nij  is  a  Markov  renewal 
process  for  which  inference  h2is  been  studied  by  Gill  (1980). 

The  results  of  the  paper  will  be  developed  for  the  nonlinear  semimartingale  regression  model 

X{t)  =  X(0)  +  r  Y{t)  a(s,  Z[a))d$  +  M[t),  (1.3) 

Jo 

where  M  is  a  square  integrable  martingale  and  Y,Z,a  are  as  before.  This  includes  diffusion 
processes  as  well  as  the  counting  processes  mentioned  above.  An  estimator  A[t,z)  of  the  time- 
integrated  conditional  “hazard”  function 

A(t,z)=  I  a{s,z)d8 

Jo 

will  be  introduced.  For  counting  processes  our  estimator  A  coincides  with  the  Nelson- Aalen  esti¬ 
mator  if  the  covariate  process  Z  is  constant.  A  kernel  estimator  d  of  or  will  be  obtained  from  A, 
as  was  done  by  Ramlau-Hansen  (1983)  for  the  Nelson- Asden  estimator. 

The  estimators  d  and  A  are  defined  in  Section  2.  Consistency  and  asymptotic  distribution 
results  for  d  and  A  are  given  in  Section  3.  In  Section  4  we  derive  confidence  bands  for  A{‘,z)  at 
any  fixed  level  z  of  the  covariate  and  introduce  a  Kolmogorov-Smirnov  type  statistic  for  testing 
the  hypothesis  that  A{-,Zi)  and  A(', Z2)  coincide  (equivalently  q:(',zi)  and  a{-,Z2)  coincide)  at 
different  levels  zi,  z^.  Technical  lemmas  used  in  the  proofs  of  the  main  results  are  given  in  Section 
5. 


2.  The  Estimators. 

{n,7,P)  will  denote  a  complete  probability  space  and  [7t,t  E  [0,1])  a  nondecreasing  right- 
continuous  family  of  sub-<T-fields  of  7  such  that  7o  contains  all  F-null  sets  in  7.  All  processes  are 
indexed  by  f  G  [0,1].  The  process  M  =  [M{t),7t)  is  assumed  to  be  a  square  integrable  martingale 
with  mean  zero  and  paths  which  are  right-continuous  on  [0,1)  with  left  limits  on  (0,1].  Suppose 
that  the  qusulratic  characteristic  (Af)  of  M  has  the  form 

{M){t)=  f%{a,Z{8),Y{s))ds,  (2.1) 

Jo 

where  7  is  a  bounded,  measurable  function.  The  covariate  processes  Y  and  Z  are  assumed  to  be 
predictable  and  V  is  an  indicator  process.  For  simplicity,  Z  is  supposed  to  be  scalar  valued.  We 
assume  that  the  processes  X,Y,Z  and  M  are  related  by  the  equation  (1.3)  which  can  be  written 
in  differential  form 

dX{t)  =  r  (t)  a{t,  Z{t))  dt  +  dM{t), 

where  a  is  a  bounded  Borel  function.  In  the  counting  process  czise 


(2.2) 


In  the  diffusion  process  C2ise  (without  censoring)  we  have  Y{t)  =  1,  Z{t)  =  X(t), 


Mit)=  r  a{a,X{s))dW[s), 
Jo 


where  is  a  Wiener  process,  (T^{t,z)  is  the  infinitesimal  variance  of  the  diffusion  and  oi{t,z)  is 
the  drift  function. 

In  order  to  define  the  estimators  a  and  A  we  need  the  following  notation.  For  z  G  R,  J, 
denotes  an  interval  of  length  ii;„  containing  z,  where  u;„  -+  0  as  n  — ►  oo.  Let  [Xi,Yi,Zi,Mi),i  = 
l,...,n  denote  n  independent  copies  of  the  generic  processes  X,Y,Z,M  which  satisfy  model  (2.1), 
(2.2).  Assume  that  Xi  and  Yi  are  observable  continuously  over  the  time  interval  [0,1]  and  Zi{t)  is 
observable  at  least  when  Yi{t)  ^  0.  Define 


=  E  r  I{Zi{s)  e  W{s)dX,is), 


(2.3) 


Y^^'>{t,z)  =  Y,nMt)€l.}Y{t).  (2.4) 

«=i 

As  an  estimator  of  A  we  propose 


1 

r('*)(s,z) 


X^’'Hds,z), 


where  1/0  =  0.  Also,  for  t  €  (0, 1)  set 


=  A{ds,z), 

where  if  is  a  bounded,  nonnegative  kernel  function  with  compact  support,  integral  1  and  6„  >  0 
is  a  bandwidth  parameter,  6„  — >  0. 

We  note  that  for  the  above  estimation  of  A  at  a  fixed  z  the  processes  X,  Z  only  need  to  be 
observed  at  times  when  Z  belongs  to  the  neighborhood  J*  of  z.  We  can  show  that  A  and  d  yield 
asymptotically  well  behaved  estimators  of  A  and  a  when  we  shrink  (i.e.  let  u;„  — *  0)  and  let 
6„  — »  0  at  appropriate  rates  as  the  sample  size  increases.  If  K  is  left  continuous  and  of  bounded 
variation,  then  by  integration  by  parts  (see  Dellacherie  and  Meyer  (1982),  Chapter  VllI,  (19.4)) 
for  n  sufficiently  large 


almost  surely  for  all  t  €  where  0  <  <  ta  <  1-  Thus,  for  n  sufficiently  large,  we  can  choose 

a  Lebesgue-Stieltjes  version  of  the  process  {a{t,z),  t  €  [ti,t2] ): 

(2-5) 

This  version  of  d  is  used  in  Theorem  2(c). 


3.  Main  Results. 

We  shall  consider  estimation  of  A{t,2),  a{t,z)  over  0<t<l,  0<z<l.  Let  C  be  a  set  in 
R  containing  U,g[o,i|  jJ"*  for  some  n  >  1.  The  following  assumptions  are  supposed  to  hold  for  all 
{t,z)  belonging  to  [0,1]  x  C. 

(Al)  For  each  t,  the  random  vector  {Z{t),Y{t))  is  absolutely  continuous  with  respect  to  the 
product  of  the  Lebesgue  and  counting  measure.  Denote  the  corresponding  density  by 
Also,  suppose  that  for  fixed  z,y  this  density  is  integrable  in  t. 

(A2)  /z(t)y(t)(z,  1)  is  bounded  away  from  zero. 

(A3)  /z(t)y(t)(z,  1)  is  continuous  as  a  function  of  t  and  z  for  each  fixed  y. 

(Bl)  a,  7  are  continuous  functions  of  t  and  z  for  each  fixed  y. 

(B2)  a  is  Lipschitz,  i.e.  there  exists  a  constant  K  such  that 

|a(ti,  zx)  -  a(t2,  Zi)  I  <  iflKti  -  t2,zi  -  zz)!! 
for  all  tx,t2,zx,Z2,  where  ||  •  [j  denotes  the  Euclidian  norm  on  R^. 


THEOREM  1.  (a)  Suppose  that  A1-A3,  Bl  hold  and  nwn  —*  oo  as  n  —*  oo.  Then 


sup  E  sup  I  A(t,  z)  -  A{t,  z)|^  — ►  0 

z  t 


as  n  —*  oo. 


(b)  Suppose,  in  addition,  that  B2  holds  and  nw^  — ►  0  as  n  — ►  oo.  Then 

nwnE\A{t,z)  -  A{t,z)'^  —*  f  h{s,z)ds 

Jo 

uniformly  over  [t,  z)  €  [0,  l]^  as  n  — *  oo  and 

lim sup  nu;„E  sup  I A(t,  z)  -  A(t,  z)|^  <  4  /  h[s,z)ds 

n— *00  t  Jo 

uniformly  over  z  6  [0, 1]  as  n  — ►  oo,  where 


(3.1) 


(3.2) 


(3.3) 


,,  7(s,i,  1) 

f  /»  iV 

fZ{,]Yl,)(z,  1) 


(3.4) 


Proof.  Define 


=  T  f  I{Zi{3)  e  h}Yi{s)dMi{s) 

.=1  -^0 

n 

=  5;/{Z.(0  e  I.}Yiit)a{t,Ziit)) 


=  f^I{Ziit)  e  I.}Yi{t)^{t,Ziit),Yiit)). 


It  follows  from  (2.1)'(2.4)  that 


X(”)  {dt,  z)  =  {t,  2)  dt  +  Af  {dt,  z) 


d{M^-\-,z))it)  =  't(-){t,z)dt. 


The  Doob-Meyer  decomposition  of  A  is 


Therefore 


where 


|A(t,2)  -  A(t,2)P  =  Ii{t)  +  hit)  +  h{t), 

hit)  =  2Iih)Iih). 


and  by  Lemmas  4  and  5 


£;SUPA(0<|  g|K|n|(,;,)  * 


"  [fSM  -  “<-)r = ”'» 


Lr'")(3,2) 

uniformly  in  s,  2  as  n  — ►  oo  if  a  is  continuous.  Also 


V 

? 

jJ: 

• 

r 

is 

? 

I 

to 


g; 


•  A  A 


(3.6) 


uniformly  in  s,  z  as  n  -+  oo  if  a  is  Lipschitz  and  nto®  — »  0.  Next,  by  Doob’s  inequality 

where  the  r.h.s.  of  the  last  equality  tends  to  4  h(s,z)ds  uniformly  in  z  eis  n  — »  oo  by  Lemma  6. 
This  proves  (3.1), (3. 3).  To  show  (3.2)  we  observe  that 

This  completes  the  proof. 

The  next  result  can  be  viewed  as  the  analogue  of  Theorem  1  for  the  estimator  a. 

THEOREM  2.  (a)  Suppose  that  A1-A3,  B1  hold  and  6n  ~  tt;„,  — ►  oo  as  n  — >  oo.  Then 

E  [Q!(t,  z)  -  a(t,  z)]^  — >  0 
for  every  t  €  (0, 1)  uniformly  in  z  as  n  — ►  oo  and 

f  E[a{t,z)  —  a{t,z)]^  dt -*  0 
Jo 

uniformly  in  z  as  n  — *  oo. 

(b)  Suppose,  in  addition,  that  B2  is  satisfied  and  nw*  — ►  0  ais  n  -+  oo.  Then 

nw^E[d{t,z)  -  a(t,z)p  -+  Kk(3,z) 
for  every  t  €  (0,1)  uniformly  in  z  as  n  — »  oo  and 

nw^  f  E[d(t,z)  -  a(t,z)]^  dt  K  f  h(t,z)dt 

Jo  Jo 

uniformly  in  z  as  n  — >  oo,  where  k  =  ^^(“) 

(c)  Suppose  A1-A3,  Bl,  B2  hold,  K  is  left  continuous,  of  bounded  variation  and  nu;„  — ♦  oo, 
nui^  — ►  0,  nwnh'^  — ►  oo.  Let  0  <  ti  <  <2  <  1-  Then,  for  the  version  of  d  given  by  (2.5), 

supE(  sup  |Q(t,z)  -  a(t,z)|)  =  C>(  ■)+0(6„). 


Proof.  From  the  decomposition  (3.5)  it  follows  that 


a(t,z)  -  a(t,z)p  =  [A(0P+  [/2(0P  +  2[/i (0/2(01, 


where 


Now 


't  —  a\\'^  7*"l(s,z) 

(r('*)(s,z)> 

Under  the  hypothesis  of  part  (a)  it  follows  from  Lemma  1  and  Lemma  3  that  nw‘f^E  is 

uniformly  bounded,  hence  E[Ii{t)\^  -*  0  and  ElIi(t)Pclt  — »  0  as  n  — ►  oo.  Under  the  hy¬ 
pothesis  of  part  (b)  it  follows  from  Lemma  8  that  nw'^E  — >  K.h{t,z)  for  all  t  G  (0,1) 

uniformly  in  z  as  n  — ►  oo.  The  bounded  convergence  theorem  implies  that  ntu^  E  [/i(0]^  dt  —* 
K  fQh(t,z)dt  uniformly  in  z  as  n  — ►  oo.  By  the  Cauchy-Schwarz  inequality,  E|/i(t)/2(0l  < 
(E  [/i(t)]^  E  1/2(0]^  the  proof  of  parts  (a),  (b)  will  be  complete  if  we  we  show  that 


F;[/2(or  = 


2  I  0(1)  under  the  conditions  of  part  (a)  for  all  t  €  (0, 1)  uniformly  in  z 


0(11;^)  under  the  conditions  of  part  (b)  for  all  t  G  (0, 1)  uniformly  in  z 


and  similarly  for  E  dt.  This  is  done  in  Lemma  7.  To  prove  (c)  we  observe  that  d(t,  z) 

a{t,  z)  =  [^{t,  z)  -I-  I^{t,  z),  where 


h{t,z]  =  Y  {A[s-,z)-  A{s-,z))dK(^-j^'^ 
^4(^,2)  =  ^^  A{s-,z)dK{^-Y^)  -  ot{t,z). 


Now 


\Iz{t,z)\<[2/hn)V{K)  sup  |A(s-, z)  -  A(s-, z)|, 

0<9<1 


where  V[K)  is  the  total  variation  of  K.  Thus,  by  Theorem  1(b) 

sup/7  sup  i/3(f,z)|  =  o(—^  ). 

*  0<<<1  yy/nWnbn' 

By  integration  by  parts,  for  n  sufficiently  large 

7 


m 


U 


.V 


'a 

S 

I 

ix*\ 


=  A((is,2) -a(t,z)  =  a(s,z)ds  -  a(i,z). 

Therefore,  by  B2,  l/4(t,z)|  =  0(6„)  uniformly  in  {t,z)  €  [ti,t2]  x  [0,1]-  This  completes  the  proof 
of  the  theorem. 

REMARK.  The  conditions  on  Wn,  bn  of  part  (c)  of  Theorem  2  are  satisfied  for  the  sequences 
bn  =  n~^  and  u;„  =  n~^,  where  1/3  <  6  <  1  and  0  <  0  <  {1  -  S)/2.  ' 

Asymptotic  distribution  results  for  A  will  be  established  under  additional  assumptions  on  X 
that  will  make  functioned  central  limit  theorems  for  martingales  easily  applicable  to  the  martingale 
part  of  (3.5).  In  the  sequel  Z)[0, 1]  denotes  the  space  of  real  valued  functions  on  [0,1]  which  are 
right  continuous  on  [0,1)  and  whose  left  limits  exist  on  (0, 1],  equipped  with  the  Skorohod  topology. 
Also,  i?[0,  ij*  denotes  the  product  space  of  t  copies  of  iP[0,  ij.  For  an  account  of  weak  convergence 
in  £>[0,1]  we  refer  to  Billingsley  (1968). 

THEOREM  3.  Suppose  that  X  has  continuous  sample  paths  or  is  a  counting  process,  A1-A3, 
B1-B2  hold  and  — ♦  oo,  nw^  — ►  0  as  n  — +  oo.  Then  for  zi,...,Z{  €  [0, 1],  all  distinct,  the  process 

(y/^(A(t,2r)  -  A(t,Zr)),te  [O,  l])'^i 
converges  weakly  in  D[0, 1]*  as  n  -+  oo  to  the  Gaussian  process 


(f/(t,z.),te[0,  l])t,,,i 


with  zero  mean  and  covariance  function 


rti  Af  3 

Cov(l7(tl,Zr,),l/(t2,Zrj))  =Sr,,r,  /  h(s,Zrj)ds 

Jo 


(where  5ri,r2  denotes  the  Kronecker  symbol). 
Proof.  By  decomposition  (3.5) 


where 


y/nWn{A{t,Zr)  -  A{t,Zr))  =  (/g (()  +  /4(t)), 


Observe  that 


nWr,E  8Up(/3(t))^  <  nw„  E  [•y(n)|^’^*^|  -  “(«.  ■^r)]  <fs  —  0 
as  n  — ♦  <30  by  (3.6).  Therefore  it  will  be  sufficient  to  show  that 

(3.7) 

weakly  in  D[0, 1]*  as  n  — »  oo.  If  X  has  continuous  sample  paths,  so  does  the  square  integrable 
martingale  I^[t)  and  by  Liptser  and  Shiryayev  (1980)  it  will  be  sufficient  to  show  that 


f  h{s,z)ds 
Jo 


for  all  t,Zr  as  n  — ♦  oo.  Since 


(/  WTOT - ^M^’''>id3,Zr)){t)  =  f  ds, 

\Jo  Yi^)is,Zr)  Jo  {Y(’-){a,Zr)f 

(3.8)  follows  from  Lemma  6.  If  X  is  a  counting  process,  by  Rebolledo  (1978)  we  will  have  to  verify, 
in  addition,  the  Lindeberg  condition 

7('‘)(s,Zr)  T(  j -  1  ■>  J  P  r. 

ntu„  /  - - ’—s  I{\/nwn,;,  w - r  >  e}ds-4  0, 

Jo  (y('‘)(s,z,))'*  '*y('»)(s,z.)  ^ 

for  all  e  >  0.  But  by  application  of  the  Cauchy-Schwarz  inequality 


<{£[ 


-  1“  L— "y('*)(s,z,)J  J  l“L  n«;„  J  ‘  '“'’*y(-)(s,zO  ' 

By  Corollary  1  and  Lemma  3  the  r.h.s.  above  tends  to  zero  for  all  Zr  if  for  all  s,  z^ 

This  follows  from  Lemma  3  and  the  Markov  inequality.  So  far  we  have  only  shown  weak  convergence 
in  Z?[0, 1]  for  each  z^.  Note  that  (/4(t))|._i  is  a  vector  of  square  integrable  martingales  that  are 
orthogonal  whenever  n  =  0,  which  is  true  for  sufficiently  large  n.  Therefore  (3.7)  follows 
from  the  previous  2irguments  by  application  of  the  Cramer- Wold  device. 

The  next  theorem  gives  the  Jisymptotic  finite  dimensional  distributions  of  the  estimator  a. 

THEOREM  4.  Suppose  that  X  hats  continuous  sample  paths  or  is  a  counting  process,  A1-A3,  Bl, 
B2  hold  and  6n  ~  ti'n,  oo,  — »  0  2is  n  — »  oo.  Then  for  all  zi,. . .  ,zi  G  [0,  . .  .  ,tk  £ 

(0, 1),  all  distinct, 


{s/nw„{a{tj,Zr)  - 

converges  in  distribution  to  the  Gaussian  random  array  mean  zero  and 


covarieince 


Cov(V’(tyj  ,  Zrj  ))^(tj2> '^'■2))  —  ^ji.J2  ^’"l >'‘2  )• 


Proof.  By  decomposition  (3.5) 


v/nu;„(a(ty,z,)  -  a{tj,z^))  =  y/nw^ih  +  h), 


where 


h=^  f  - -M^^\ds,Zr). 

bn  Jo  ^  bn  ^Yi’')a,Zr) 

It  follows  from  Lemma  7  that  y/nwnIs-^0  as  n  —*  oo.  Therefore  it  will  be  sufficient  to  show  that 


Now  /e  =  /6(l),  where 


hi 

and  (f6(»'))jli,r=i  is  array  of  kl  squeure  integrable  martingales  that  are  orthogonal  for  n  suffi¬ 
ciently  large.  If  X  has  continuous  sample  paths  it  follows  from  Remark  2  in  Liptser  and  Shiryayev 
(1980)  suid  the  Cramer- Wold  device  that  we  only  need  to  show  ni(;„^(/6)(l)-^  K.h[tj,Zr)  as  n  -+  oo 
for  all  tj,  Zp.  This  is  done  in  Lemma  8.  If  X  is  a  counting  process  we  have  to  verify,  in  addition, 
that 

Jo  ^bn  \  bn  J'  (y('‘)(5,z))^  '  bn  ^  /y(")(s,z)  i 

Eis  n  — ►  OO  for  all  e  >  0.  As  in  the  proof  of  Theorem  3  this  follows  from 

P[V^w„  i  K{i^)  >  e|  <  >  e| 

(for  some  e  >  0,  since  K  is  bounded) 


This  completes  the  proof  of  the  theorem. 


(by  Lemma  3  and  the  Markov  inequality). 


4.  Confidence  Bands  and  Hypothesis  Tests. 

In  order  to  use  the  previous  theorems  for  inference,  an  estimate  of  h{t,  z)  (as  defined  in  (3.4)) 
and  of  /q  h{3,z)ds  =  H{t,z)  is  needed.  The  following  theorem  provides  consistent  estimators  for 
both  of  these  quantities  in  the  case  that  X  is  a  counting  process. 


THEOREM  5.  Suppose  that  A1-A3,  Bl,  B2  hold,  A”  is  a  counting  process.  Define 


H{t,z) 


1 


I 


A<'*)(ds,z), 


for  t  G  [0, 1], 

hit,z)  =  ^  H[d3,z), 

for  t  G  (0, 1).  If  bn  ~  nWn  -*  oo  as  n  cx  then 


E\h{t,z)  -  h{t,z)\  -*  0 

for  all  t  G  (0, 1)  uniformly  in  z  €  [0, 1]  as  n  — ♦  oo.  If  nu;„  -+  oo  then 


E  sup|.H(t,z)  -  H{t,z)\-*0 

t 

uniformly  in  z  G  [0,  l]  as  n  — ►  oo. 

Proof.  By  decomposition  (3.5) 


\h{t,z)  -  h{t,z)\  <  hit)  +  hit), 

where 

. . . I 

'■"’■Itt/.' 

By  Lemma  8  (with  replaced  hy  k  K)  and  ~  we  have 


r" 


PJTWWVJlKnK'V  vrsi'y  'k-v-.v.-w  v-«  v’v'.v^-j-.v-rsrj-rvs^-v-.'srL^.'v-.'r.-r.'v'- ».-.  v  %-  ’^.-. ».-  •-  > 


SS 

N’ 

lV 


:*S 


•y 

:? 

V 


I 

?■ 

I 

s 


[/2W]'^  <-;K  [  ( V^)  ^  da  ->  0 

^  b^^Jo  ^  hr,  J  {YM{a,z))* 


for  all  t  e  (0, 1)  uniformly  in  z  2a  n  — ►  oo  by  Lemmas  1,3.  Similarly 


\H{t,z)  -  H{t,z)\=  nw„  f  ^ — ^X<")(da,z)-  /  /i(s,  z)  da  <h{t)  +  h{t), 

Jo  (^''‘/(a.z))  Jo 


where 


h{t)  =  nw^  f  da  -  f  h{s,  z)  da  I , 

Jo  iYM{s,z)f  Jo  ^  ^  I 

h{t)=  nwn  f  — ; — ^ - r-M^"^(da,z) 

Jo  (y('‘)(a,z)f  ^  ^ 


But 


sup  sup /a  (t)  <  /  3up£? 
X  t  Jo  x,t 


nw„ 


a*'*){a,z) 

(yl^lCa.z))' 


-  h{s,z) 


ds  —*  0 


by  Lemma  6  (with  7^"^  replaced  by  a*"))  and 

sup£:sup[/4(t)]^  <  f  supJS  [(nwn)^ 

*  ‘  Jo  ».*  I-  (y'’*)(s,z)) 

as  n  — >  00  by  Doob’s  inequality  and  Lemmas  1,  3. 


da  — ►  0 


In  the  diffusion  process  case,  in  which  <T^{t,z)  is  assumed  to  be  known,  the  following  theorem 
provides  consistent  estimators  for  h{t,z)  and  H{t,z). 


THEOREM  6.  Suppose  A1-A3,  Bl,  B2  hold,  X  is  a  diffusion  process.  Define 

[L 


for  t  G  [0, 1], 


for  t  €  (0, 1).  If  6„  ~  u;„,  nw^  — ♦  00  then 

E|h(t,z)-h(t,z)|  -0 

for  all  t  e  (0, 1)  uniformly  in  z  G  [0, 1]  as  n  — ►  00.  If  nWn  —*  00  then 
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EG6S4rfl&0‘j'2^.;)^>  .vC-<.y2V,L-'flA  / 


uniformly  in  z  €  [0, 1]  as  n  — ►  oo. 


E  sup  |^{t,z)  -  H{t,  z)\  —y  0 

t 


Proof.  It  follows  from  6„  ~  ti;„  and  Lemma  8  (where  we  replace  by  [s ,  z)  Y {s ,  z) 

and  hy  k  K)  that 


E\h{t,z)  -  h{t,z)\ 


yt'*)(s,z) 


ds  —  h{t,  z) 


—  0 


for  all  t  G  (0, 1)  uniformly  in  z.  Next 


E  sup|.ff(f,  z)  —  H{t,z)\  <  sup  E 

t  »yZ 


nw„ 


^is,z) 


yi^iCs.z) 


his,z)  , 


which  tends  to  zero  by  Lemma  6  with  ^^"^(sjz)  replaced  by  a^(s,  z)  z).  This  completes  the 

proof  of  the  theorem. 


Confidence  bands  for  A(-,z). 

Under  the  conditions  of  Theorem  3 


H{^z)  +  z)  ~*^°(h(1,z)^+  ff(t,  z) ) 


weakly  in  Z)[0, 1]  as  n  — ♦  oo,  where  is  the  Brownian  bridge  process.  By  application  of  Theorem 
5  in  the  counting  process  case  or  Theorem  6  in  the  diffusion  process  C2ise  we  obteun  the  following 
eisymptotic  100(1  -  a)%  confidence  band  for  A(-,z): 


where 


A(t,  z)  ±  Co 


1  ,  H{t,z) 
H{\,z) 


),  te[o,i], 


P[  sup  (H^°(t)|  >  Co]  =  a. 

«€[0,l/2l 

A  table  for  the  distribution  of  sup<g|o,i/2|  |^”(0I  found  in  Hall  and  Wellner  (1980). 

Testing  equality  of  A  at  two  different  levels  of  the  covariate. 

We  now  introduce  a  test  statistic  for  testing  the  null  hypothesis  Hq  ■  zi]  =  A{t,  Zo)  for  all 
t  €  [0, 1],  where  Zi,  Z2  we  two  prechosen  values  of  z.  Define 


Ai2(«)  =  A(t,zi)  -  A(t,Z2) 
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Ai2it)  =  A{t,zi)  -  A{t,Z2) 
Hi2{t)  =  Hit,zi)-  H{t,Z2) 
Hx2{t)  =  Hit,zi)-H{t,Z2). 
Then  under  the  conditions  of  Theorem  3  we  have  that 


weakly  in  i?[0, 1]  m  n  — »  oo.  Set 


(^12(0  -  -AiaCO  ) 


WO 


(  Hi2(t)  ^ 
'^12(1)  +  Hi2{t)' 


I  =  \/ntu„^i2(l) 


sup 

t€|0,l| 


-^12(0  ~  -^12(0 
^12(1) +^12(0 


Then  in  the  counting  process  and  diffusion  process  cases  considered  above  as  n  — >  00,  where 
sup4g[o_i/2|  Therefore  an  asymptotic  test  of  size  ot  can  be  carried  out  by  rejecting  Ho 

if  and  only  if  |  >  c^,  where  P(^  >  c,,)  =  a.  Finally  we  mention  that  Theorem  4  can  be  used  to 
construct  asymptotic  tests  as  in  Rao  (1973)  for  testing  equality  of  a  at  any  finite  number  of 
vzilues  of  t  and  z. 


5.  Technical  Lemmas. 

LEMMA  1.  Suppose  that  Al,  A3,  B1  hold  and  n«;„  — ►  00  as  n  — »  00.  Then 


nw„  J 


for  all  nonnegative  integers  k,  uniformly  in  s,  z  as  n  — ♦  00,  where 


g(s,z)  =  /z(,)y(,)(z,  1)7(3, z,l). 


(5.1) 


Proof.  By  the  multinomial  theorem 

Since  |^»(s)  -  z|  <  «;„  implies  |7(s,  Zi(s),  1)  —  7(3, z,l)|  <  e„  uniformly  in  s,  z  for  some  e„  — ►  0  as 
n  — >  00,  we  have  for  ji  ^  0 


Eiiiz,(s)  e  J.}y;(s)7(3,-^.(s),y.(s))f  =  £;[/{z.(s)  e  j,}y;(s)(7(3,z,i)  +  o(en)]'’ 
=  ( ( (7(s,  1)  y-  +  0(c„) )  E  [/{Z. (s)  e  7,}  V,(s)  j  =  /, . 
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Also,  by  uniform  continuity  of  /, 


\E[I{Zi{s)  e  /,}yi(s)]  -  runfz{,)Y[,){zA)\  <  f  |/z(.)y(.)(u,  1)  -  fz(,)Y(,){z,l)\du 

J  I, 


uniformly  in  s  and  z.  Therefore 


h  =  {'l{s,z,l)y'Wnfz{,)Y(»){z,  l)  +  iu„o(l) 


uniformly  in  s,  z  and 


nil’*’]  " 


=  )*  +  o(l). 


COROLLARY  1.  Suppose  that  Al,  A3  hold  and  nu;n  — »  oo  as  n  — +  oo.  Then 


®  1“;^ '>1 


as  n  — *  oo  for  all  s,  z,  it  >  1, 


uniformly  in  s,z  as  n  —*  oo  and 


uniformly  in  s,z  as  n  — ►  oo. 


—  0 


■  z 

—  0 


LEMMA  2.  Suppose  X  ~  binomial  (n,p),  0  <  p  <  1.  Let 


Then  for  each  integer  ^  >  1 


1/X,  if  X  >  0; 
0,  ifX  =  0. 


\  nv  / 


r 


WVVV'.'VVWr.V.».''.V  vr^'.-v;  vt-.,-  •rj^.’srjrawVTyrni-Tjr-jr^ -U>  V irw  <r»  « - «.- . -.-  -u 


Proof. 


'•1*  =  V'  1  - ^ _  n‘  ««-  =  V  1  («  +  l)---(«  +  fc)  _J±_  .  „„-. 


(■  +  jfc)!  (n-t)!^*' 


-Z^T^TiT^rroiP’  - (T 

1=1  '  '  i=l  ' 

<-p^  '’■  ’"-■  =  ^  t 


LEMMA  3.  Suppose  Al,  A2  hold  and  nti;„  — ♦  oo  as  n  — »  oo.  Then 


uniformly  in  s,z  as  n 


.r(")(s,2)J 

oo  for  every  integer  Jfc  >  1,  where  1/0  =  0. 


Proof.  Set  m  =  inf,^,  /2(«)y(») (■z.  !)•  Then  y'^'*^(s,2)  has  a  binomial  distribution  with  parameters 
(n,pl”l(s,z)),  where  p*")(a,z)  >  mwn,  ao  the  previous  lemma  applies. Therefore 


nwn  1*  ^  /(A:  +  l)nu>n\fc 
y('*l(s,2)i  “  V  nmwn  ' 


In  the  following  lemma  we  will  use  the  notation  J*'*^(s,z)  =  /{y*'*)(s,2)  ^  0}. 


LEMMA  4.  Suppose  that  Al,  A2  hold  and  nw„  — *  oo  as  n  — ►  oo.  Then 

Ejl  -  J<")(s,2)l*'  <  exp{-nu;„  inf /z(,)y  (,)(z,  1)} 

«,Z 

for  each  integer  A:  >  1. 

Proof.  Set  m  =  inf,,*  /z(,)y(,)(z,  1).  Then  m  >  0  and 

Ell  -  =  p («/<"’ (a, z)  =  0)  =  (1  -  P[Z[s)  e  7*.  y(s)  =  1])" 

<  (1  -  miUn)”  <  exp{-mnui„}. 


LEMMA  5. 


a*"! (s,  z) 

yi'*)(s,2) 


a(s. 


o(l)  uniformly  in  s,2  if  a  is  continuous 

0{w„)  uniformly  in  s,2  if  a  is  Lipschitz. 


Proof.  By  definition  of  a^")(s,2) 

where  =  o(l)  (=  0(u;„) )  uniformly  is  a,  z  if  a  is  continuous  (if  a  is  Lipschitz). 

LEMMA  6.  Suppose  that  A1  -  A3,  B1  hold  and  nti;„  -+  oo  as  n  — *  oo.  Then 


sup  E 


nwti 


where 


h{3,z)  = 


(y('‘)(s,z)) 


-  Hs,z) 


0, 


Proof. 


where 


nw„ 


(y(")(s,2)V 


-  h{s,  z) 


<  h  +  I2, 


By  Lemma  4,  /o  — *  0  uniformly  in  s,z  as  n  —♦  00.  Now  by  application  of  the  Cauchy-Schwarz 
inequality 


where 


h  =  E 


7'")(s,z)  /r<'*)(s,z)\2 


nw„ 


V  nw„  / 


Also  4/3  <  A  +  /&,  where 


U^E 


nw„ 


g{s,z) 


h  =  C^E 


/y'")(s,z)N^2 

I - )  - 

\  nWn  ' 


where  C  =  sup,  ^  h(3,z).  But 
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ii.  I.  ii.t 

•  M 


But  sup,  ,  E{1  —  z))^  =  0(l/(nwn)  )*“  for  all  integers  fc  >  1  by  Lemma  4.  This  and  6„  ~  tti„ 

imply  Wn~^l2{t)  =  0{l/{nwn^^^’‘)  )*'  for  3/k  <  6  uniformly  in  t,z  as  n  — ►  oo.  Thus  l2{t)  =  0{w'^) 
and  Jq  hit)  dt  =  0(u;^),  Finally 


hit)  = 


o(l)  for  all  t  €  (0, 1)  uniformly  in  r  if  ot  is  coninuous, 


[  Oiw^)  for  all  t  €  (0, 1)  uniformly  in  z  if  a  is  Lipschitz. 
Uniform  boundedness  of  hit)  and  the  dominated  convergence  theorem  imply 


r  hit)  dt= 

Jo 


o(l)  uniformly  in  z  if  a  is  continuous, 

Oiw^)  uniformly  in  z  if  a  is  Lipschitz. 


This  proves  the  lemma. 


LEMMA  8.  Suppose  that  A1-A3,  B1  hold  and  nw^  — »  oo  as  n  -+  oo.  Then  for  each  r  6  [0, 1] 

sup  (^)  -  0 

as  n  — *  oo  for  all  t  €  (0, 1),  where 

0,  if  r  <  t; 

ifr  =  t; 

.  K,  if  r  >  t. 

Proof.  For  t  <  t  the  theorem  is  obvious.  Suppose  t  >  t.  Then 


^  "io  V  6n  ))  (r(")(s,z))3 


ds  -  KiT,t)hit,z)  <  hit)  +  hit), 


where 


^2(0=  ^  /i(s,z)cis  -  «(r,t)/i(t,s)  . 

Note  that  ui„/6„  — »  1.  Therefore  hit)  — *  0  for  all  t  €  (0, 1)  uniformly  in  z  as  n  — ►  00  by  Lemma  6 
and  hit)  — ►  0  for  all  t  €  (0, 1)  uniformly  in  z  as  n  — ►  00  by  continuity  of  h. 


-  his, z)  ds 
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